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Fluid studies PVT -Summary 


► Preface 

► Reminders on thermodynamics 

• Petroleum genesis 

• Pure component and mixture equilibrium 

• Hydrocarbon fluid classification 

• PVT experiments 

• Constant volume gas depletion 

• PVT Modeling 

► Thermodynamic model 

• Liquid-Vapor equilibrium 

• Equation of state (EOS) 

• Compositional grading 

• Matching of experimental data 

► Characterization and modeling of heavy fractions 

• Heavy fractions 

• Viscosity modeling 



► Measurements 

• Sampling 

• Analysis 

► Fluid synthesis 

► Gas injection 

• Specific PVT experiments 

• Miscibility (FCMP, MCMP, Kr...) 

• Compositional matching 

► Interfacial Tension 

► Calculation of OOIP/OGIP 

• Black-oil data 

• Compositional data 

► Fluid synthesis examples 

► Fluid studies quick sum up 

IFP Training I 4 

Fluid studies PVT ^ 







Preface 




Sampling Miscible/lmmiscible 

displacement 
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We need to know 


1. The composition of the production well stream and its temporal variation 

2. The separator specifications including LPG 

3. The design of the completions by identifying the spatial fluid distributions in 
the vicinity of the wellbore 

4. Gas injection/re-injection 

• Identifying gas composition 

• The interaction of the injected gas and the reservoir fluid 



AND 


5. The ultimate recoveries of components, under different drives, mixing/no 
mixing, single depletion, etc. 

6. The amounts and composition of liquids left in the reservoir not recovered 
(especially in gas condensate reservoirs) and their properties: density, 
Surface Tension, viscosity. 

7. To detect spatial variation of the PVT properties 

8. To identify and adjust data inconsistencies 
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AND 


For Reservoir Performance simulation software: 

9. To create "Black-Oil" PVT tables for Black Oil Simulation 

10. To create Compositional PVT models for compositional simulation 



OBJECTIVE OF GATHERING AND ANALYZING PVT DATA 


► To calibrate models (Equations of State, EOS) 

► Why? 

• We cannot measure all the characteristics of hydrocarbon fluids. EOS provides one 
consistent source of PVT data 

• Experiment problems (cost, reliability, accuracy and precision) 
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PVT 


► Sometimes people spend a lot of time on the rock characteristics and put less 
emphasis on fluid PVT. 

► Fluid PVT is as important as rock properties as it is directly related to reserves 
and the dynamics of the reservoirs for performance estimations 
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ROCK 
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Petroleum genesis 


2. Reminders on 
thermodynamics 




Petroleum genesis 




TRANSPORTATION: T=ld 




PRODUCTION 

T=1Y 


RESERVOIR 

T=1000Y 


STORAGE 

T=1Y 


Organic Material 



BASIN FORMATION 
T=100MY 


SOURCE ROCK 
CATAGENESIS 
T=1MY 


Kerogen 



Fluid Studies - PVT 


!PP Training 


15 


Petroleum genesis: Hydrocarbons 




DR( CARBONS 



Hydrogen: H 



Carbon: C 


O i 


Valency = 1 


2 



4 

Valency = 4 




Simple bond 




Different types of bonds 
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Petroleum genesis: Hydrocarbon classification into groups 




CARBON atom of Valency 4 HYDROCARBONS of different 

r , , I 1 size a a d shape. 

Different types of bonds ^ 



Petroleum genesis: Hydrocarbon classification into groups 
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Petroleum genesis: Hydrocarbon classification into groups 




SATURATED HC 


Normal 

Paraffins 

Isoparaffins 

Naphtenes 

Straight chains 

4 s 

k 

4$4 

4^4y 

4$4$4 > 

4f4$4y 

4$4f4f4> 

Ramified chains 

4^4 

4 ° 

4* *4^4 

4 ° 

4 ° 

0&5 

4#4^4 

4 ° 

Cycles 

CUD <2 

•U 
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UNSATURATED HC 


Aromatics 

Olefins 

Benzene cycles 

Double and/or 
triple bonds 

ix 


y 

Jyl 

V ; S= 



Ak 


Vr 1 

k 
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NORMAL and ISO ALKANES: Cn 2n+2 

► Group of hydrocarbons consisting of linear molecules with the formula 
C n H 2n+ 2 - Methane, CH 4/ is the simplest member. Higher members, starting at 
about C 18 , are wax-like and are called paraffins. 

CYCLOALKANES: Cn HLn 

► At low carbon numbers Cycloalkanes are less stable than normal alkane 
counterparts 

• thus present in much smaller amounts 

• but can have significant effect on the phase behavior 
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Petroleum genesis 




Normal boiling point of paraffins 




Iso-paraffins are lighter than 
Normal-paraffins 
Tb(l-C„H 2n J < T b (N-C„H 2n J 
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Petroleum genesis 



Aromatics 


► Aromatic hydrocarbons contain one or more benzene rings 

► Stable rings with 3 carbon-carbon double bonds - C 6 H 6 

► They have some significant effect on the phase behavior. 
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Non Hydrocarbons: Sulfur components 



Hydrogen 

sulfide 


o 


NBP= - 60°C 


h 2 s 



Gaseous, extremely 
TOXIC, Corrosive 


Mercaptans 


► Methyl-mercaptan 



NBP= 6°C 


ch 3 sh 

• Ethyl-mercaptan 



NBP= 35°C 


c 2 h 5 sh 


Acidic components, 
bad smell 


Sulfides, 

Disulfides,. 


* Di-methyl-sulfide 



NBP= 37°C 


CH3SCH3 


• Di-methyl-disulfide 





NBP= 109.6°C 
CH3S2CH3 

Non acidic, BUT 
decompose due to 
heating into mercaptans 


Thiophenic 

compounds 

• Thiophene 



NBP= 84° C 


c 4 h 4 s 


• Benzo-thiophene 



c 8 h 6 s 


Non acidic 
components 
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Non Hydrocarbons: Other components 



Non hydrocarbons 


► Nitrogen - Pyridine - Most others are unstable 


C 6 H 5 N 

* Pyridine 



► Oxygen - Aromatic alcohols, aldehydes, carboxylic - pheno1 
acids 


c 6 h 5 oh 



► Organo- Metallic -Trace amounts 
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Some definitions 


► Intensive properties 

• They do not depend on the quantity of material (density, viscosity, permeability) 

► Extensive properties 

• They depend on the quantity of material (mass, volume) 

Extensive properties can be summed but intensive properties require specific averaging 



Some definitions 

► Isothermal 

• An isothermal process takes place at constant temperature 

► Isobaric 

• Isobaric processes occur under conditions of constant pressure 

► Isochoric 

• Describes a process or system change which takes place at constant volume 

► Isenthalpic 

• Describes a process or system change which takes place at constant enthalpy 
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Pure component equilibrium 
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Vapor pressure curve of C2H6 




-40 -20 0 20 40 60 

Temperature (°C) 
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Vapor pressure curves of several normal alkanes 




Temperature °C 
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Acentric factor w 



1/Reduced temperature (1/Tr) 


1 1.25 1/0.7 1.5 



(Pitzer's factor) 
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Properties of pure components 



Name 

Mw 

(g/gmol) 

Tc 

(K) 

Pc 

(bar) 

Acentric 

factor 

methane 

16.04 

190.4 

46.0 

0.011 

ethane 

30.07 

305.4 

48.8 

0.099 

propane 

44.10 

369.8 

42.5 

0.153 

n-butane 

58.12 

425.2 

38.0 

0.199 

i-butane 

58.12 

408.2 

36.5 

0.183 

n-pentane 

72.15 

469.7 

33.7 

0.251 

i- pentane 

72.15 

460.4 

33.9 

0.227 

n-hexane 

86.18 

507.5 

30.1 

0.299 

benzene 

78.11 

562.2 

48.9 

0.212 

cyclohexane 

84.16 

553.5 

40.7 

0.212 

n-heptane 

100.20 

540.3 

27.4 

0.349 

n-octane 

114.23 

568.8 

24.9 

0.398 

n-nonane 

128.26 

594.6 

22.9 

0.445 

n-decane 

142.29 

617.7 

21.2 

0.489 

nitrogen 

28.01 

126.2 

33.9 

0.039 

carb dioxoyde 

44.01 

304.1 

73.8 

0.239 

hydr sulfide 

34.08 

373.2 

89.4 

0.081 
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Fluid Studies - PVT 


Temperature 


Mixture 


Tc 



Phase envelope 




Critical point 
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Phase envelope 




Liquid like 
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Phase envelope 
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C 2 H 6 - C 6 H 14 binary mixtures 



(U 

Q) 

±- 

3 

to 

</) 



Temperature (°C) 
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Under saturated oil 



► Point 1 indicates that the initial reservoir 
pressure is undersaturated. 


► 


► 


► 

► 


As reservoir depletion proceeds from 1 to 2, 
reservoir hydrocarbon fluid moves in the pore 
space in the form of single phase oil and gas-oil 
ratio (the ratio of gas produced at surface 
conditions to the oil produced at surface 
conditions) at surface conditions is constant. 
The amount of oil production necessary to 
bring the pressure down from point 1 to 2 is a 
small fraction of the total production. 

As production continues and pressure is 
reduced along path 2 to 3, gas comes out of 
solution in the reservoir. 

At point 3, the effective permeability to oil and 
gas are about equal. 

But the viscosity of the oil is about two-orders 
of magnitude larger than viscosity of the gas. 


Pressure path 



► Therefore, at point 3, usually more than 90% of 
the reservoir flow stream volume is gas. 


Fluid Studies - PVT 


!f=P Training 


Volatile oil 




Volatile oil 



► Volatile oils have fewer heavy molecules than 
black oils. 

► The critical temperature of volatile oil is much 
lower than that of black oil, and is close to 
reservoir temperature. 

► The figure on the left shows that, as reservoir 
pressure is reduced below bubble point 
pressure, large volumes of gas leave the 
solution. 

► This rapidly reduces the effective permeability 
to oil and the reservoir flow stream becomes 
mostly gas within a few hundred psi below the 
bubble point. 

► The effective permeability to oil may become 
virtually zero and the flow stream may 
essentially be gas long before reservoir 
reaches the point 3. 



Fluid Studies - PVT 


!PP Training 


43 


Under saturated oil 




Pressure 

“ % > 

• • 

• • 

• • 

• • 


•.Formation pressure 


Depth 


Saturation pressure 
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Saturated oil 



Saturated oils 
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Pressure 
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Condensate gas 




Condensate gas 



► Retrograde gases have even fewer heavy 
molecules than volatile oils. 


► The critical point shifts to the left and 
downward in the phase diagram and the critical 
temperature is usually less than reservoir 
temperature. 

► Retrograde condensate appears in the reservoir 
pore spaces at pressure below the dew point 
pressure. Throughout most of the reservoir, 
since the amount of liquid in the pore space is 
usually smaller than critical oil saturation the 
effective permeability to this condensate is 
zero and little is produced. 

► Along line 2 to 3, the condensate builds up first 
and then revaporizes at the lower pressures. 

► This behavior is typical for constant 
composition expansion type application. 

► At reservoir conditions can we see re- 
evaporation? 
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Condensate gas 



Pressure 

> 


Formation pressure 


Saturation pressure 


Depth 


V 
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Wet gas 



Pressure path 
in reservoir 

1 



Wet gas 


IFP Training 50 


Wet gas 



The composition of a wet gas 
contains much fewer heavy 
molecules. 

Because the phase diagram covers 
a much smaller temperature range, 
the pressure depletion path in the 
reservoir does not enter the two- 
phase region. 

The reservoir fluid is gas 
throughout the life of the reservoir. 

However, separator conditions lie 
within the two-phase envelope, 
indicating that some liquid will 
condense at the surface 



Pressure path 
in reservoir 

1 


Wet gas 
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Wet gas 



Pressure 


Depth 


‘•Formation pressure 
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Dry gas 


Pressure path 
in reservoir 

1 



Dry gas 


2 • 


Dry gas 


► Dry gas is virtually pure methane. 

► The two-phase envelope is small and 
lies below reservoir conditions and 
to the left of surface condition 


► The fluid is theoretically gas both in 
the reservoir and at the surface. 



Pressure path 
in reservoir 

1 
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Dry gas 


54 
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Dry gas 




Pressure 

> 


'•Formation pressure 


Depth 
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Fluid classification 



CF: Critical Fluid 



GC: Condensate Gas 
WG: Wet Gas 
DG: Dry Gas 


Reservoir 

conditions 


LO: Light Oil 

BO: Black Oil 

HO: Heavy Oil 
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Fluid classification: saturated or not 


Temperature 


Fluid Studies - PVT 


!f=P Training 


Fluid Studies - PVT 


Reservoir 
" ^ conditions 


Temperature 


Saturated 


Reservoir 

conditions 


Fluid classification: oils 



o> 


</> 

o 


CL 
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Fluid classification 
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Separator 

Density 

Composition 

Well test 

GOR 

of STO 

Cl 

C2 

C3 

C4 

C5 

C6+ 


(Sm 3 /m 3 ) 

(kg/m 3 ) 

(mole %) 










Heavy oil 

<10 

>900 


3 

4 

5 

8 

80 

Standard oil 

< 500 

800-900 

45 

4 

4 

3 

2 

42 

Critical fluid 

# 700 

750 -850 

55 

10 

8 

5 

6 

16 

Condensate gas 

700-800 

700-800 

75 

8 

5 

2 

2 

8 

Wet gas 

>15000 

700-800 

90 

5 

3 

1 

1 

1 

Dry gas 

infinite 


95 

3 

1 

1 














Separator 

Density 

Composition 

Well test 

GOR 

of STO 

Cl 

C2 

C3 

C4 

C5 

C6+ 


(Sft r 3 /barrel) 

(lb/ft 3 ) 

(mole %) 










Heavy oil 

<60 

>56 


3 

4 

5 

8 

80 

Standard oil 

< 3000 

50-56 

45 

4 

4 

3 

2 

42 

Critical fluid 

# 4000 

47-53 

55 

10 

8 

5 

6 

16 

Condensate gas 

4000-4500 

47-53 

75 

8 

5 

2 

2 

8 

Wet gas 

>85000 

44-50 

90 

5 

3 

1 

1 

1 

Dry gas 

infinite 


95 

3 

1 

1 
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Bottomhole properties 




GOR 

Bottomhole 

Bottomhole 

Well test 

separator 

density 

Viscosity (**) 


(Sm 3 /m 3 ) 

(kg/m 3 ) 

(cPo) 

Bitumen 


>950 

>1000 

Heavy oil 

< 10 

850 - 950 

sevl 10 2 to sevl 10 3 

Standard oil 

< 500 

550 - 850 

0.2 to several 

Critical fluid 

# 700 

450 - 550 

0.2 to several 

Gas condensate 

700 - 800 

300 - 450 

<0.3 

Wet gas 

> 15000 

100 - 400 (*) 

<0.3 

Dry gas 

infinite 

100 - 400 (*) 

<0.3 


*: highly pressure-dependent 
**: highly temperature-dependent 





Y m \ 

Evolution of density 

\md viscosity \ 

Single-phase 
(P>P sat) 

p and^p \jary(jA/ith / 

Two-phase 
(P<P sat) 

Oil 

p and p var]^contrary toj^ / J 


Gas 


p and p vary with P 
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Oil and gas specific gravity 



Stock Tank Oil Specific Gravity 

► API definition: ° API = (141.5 / SG) - 131.5 

SG being the specific gravity of oil w.r.t water @ 60°F 

• Condensates, very light oil SG < 0.8 (> 45°API) 

• Light oil 0.80 < SG < 0.86 (33-45°API) 

• Standard oil 0.86 < SG < 0.92 (22-33°API) 

• Heavy oil 0.92 < SG < 1.00 (< 22°API) 

Gas Specific Gravity 

► Gas SG is defined w.r.t air @ Std conditions (60°F, 1 atm). 

► Gas SG can be determined from the gas composition (Mw): 

SG = (gas Mw) / (air Mw) with air Mw = 28.996 g/gmol 

► SG = Mw / 29 
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PVT experiments 



► Oil: 

• Bubble point 

• Flash test 

• Multi stage separation 

• Differential liberation 

• Constant Mass Expansion (CME) 

• Viscosity under pressure 

► Gas: 

• Dew point 

• Multi stage separation 

• Constant Mass Expansion (CME) 

• Constant Volume Depletion (CVD) 
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PVT experiment: 

constant mass expansion for an oil 

Graphical determination of the 
saturation pressure 




PI 


P2 


P3 


Psat 


P4 

P5 

P6 


Volume 
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PVT experiment: flash test for an oil 




► Objective: 

• To determine reservoir fluid composition 

► Main properties derived from this experiment: 

• Molecular composition of reservoir fluid 

• Rough GOR value 

• Rough Bo value 

• Density and viscosity of oil at atmospheric conditions 



PI > P sat 
T reservoir 



v std 

gas 


V st0 

w oil 


15°C - 1 atm 


Bo = 


res 


V r< 

oil 


V 


Sto 


oil 


GOR = 


V 


Std 


gas 


V 


Sto 


oil 
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PVT experiment: flash test for an oil 



Characterization of oil and gas phases 


Gas 

► Composition (GC) 

► %H 2 S 

► Trace components 

► SG 


Oil 

► Composition (GC) 

► Heavy fraction (GPC) 

► Wax 

► WAT 

Wax Appearance Temperature 

► Pour point 

Lowest Flowing Temperature 

► Density 

► Viscosity 
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PVT experiment: multi stage test for an oil 



► Objective: 

• Simulation of field test or process separation scheme 


► Main properties derived from this experiment: 

• GOR 

• Process Rsi 


Process Boi 

Gas composition at each step 



Gas treatment 
Fuel gas 
Injection 


V 


Std 

gas 


Stock \j 
tank 


Sto 
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PVT experiment: multi stage test for an oil 



Characterization of oil and gas phases 


Gas 

► Composition (GC) 

► %H 2 S 

► Trace components 

► SG 


Oil 

► Composition (GC) 

► Heavy fraction (GPC) 

► Wax 

► WAT 

► Pour point 

► Density 

► Viscosity 
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► Objective: 

• Simulation of reservoir depletion 


► Main properties derived from this experiment: 


• Liberation GOR 

• Differential Rsb 

• Differential Bob 

• Gas composition at each step 

• Residual oil composition and 

properties 
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Reservoir temperature 


Pi<P f 


r 


T = 15 °C 



P 2 <Pi 




1 atm P = 1 atm 
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Differential vaporization 





V' 


gas 


V 


Std 


gas 



Bo = 


V 


oil 


V 


Sto 


oil 


GOR 


Total 


v std 

gas _ 

YT “ 
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Differential vaporization 





Reservoir Pressure 

> 

Bubble point pressure 
No compositional change 
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Differential vaporization versus multi-stage test 



► Differential liberation. The final volume of liquid phase remaining in the cell 
at standard conditions is called RESIDUAL OIL 

► Multi-stage test. Oil that results when one m3 or barrel of oil is flashed 
through a certain surface separator is called STOCK TANK OIL. 

► RESIDUAL OIL and STOCK TANK OIL are both products of the original oil but 
are developed by different pressure-temperature routes. 

• Multiple series of flashes at the elevated temperature of the reservoir 

• Some stage flash at low temperature and pressure 


B 


Dif 

0 




>P 


Sep 

STO 



> R 


Sep 

S 
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Viscosity measurement 



Constant Temperature = reservoir T 
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Measurement of pressure drop across a capillary tube: 


kAP 




Oil unit performances 



► Range of work: 

• Pressure: 1 bar - 1000 bar 

• Temperature: -20°C to 180°C 

• GOR: 5 to 700 m 3 /m 3 

• Viscosity: 0.2 to 200 cPo 

► Accuracy: 

• Volumes: 

• resolution 5 10" 6 cm 3 

• accuracy 10" 4 cm 3 

• Temperature: ± 1°C 

• Pressure: 

• HP: ±1 bar 

• BP: ±10' 3 bar 

• Viscosity: ±10' 2 cPo 
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PVT experiment: constant mass expansion 
for a condensate gas 

► Objective: 

• Fluid behavior at reservoir conditions 



► Main properties derived from this experiment: 

• Liquid drop out curve 

• Compressibility factor curve 

► Remarks: 

• The measurement accuracy of a very small volume of liquid is low. 

• The dew point may never be known better than +/-1.5 bars (20 psi) 

► Notice that CME = CCE 
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PVT experiment: constant mass expansion 


for a condensate gas 
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PI > Psat P2 = Psat P3 < Psat 
Retrograde 



P4 P5 P6 = Psat P7 

Normal 



PVT experiment: constant mass expansion 
for a condensate gas 



Phase envelop - Liquid deposit 



Temperature T res . Volume 
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PVT experiment: constant mass expansion for a 


condensate gas 


Compressibility factor: Z 

z=^V- 



Compressibility factor Z 




Pressure 


Optical Detection 


n 


Refracted light 


Constant T= Reservoir Temp 


pv 

Mw 


Process 



Optical density @ 202.0°C 

40 

35 

30 

25 

20 

35 

5 



















( 

) 200 

400 600 800 100C 

Pressure 



Material balance ^ p @ bhc 
Analysis ^ Mw 
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PVT experiment: constant volume depletion 
for a condensate gas 



► Objective: 

• Simulation of reservoir production for condensate gas and volatile oil 


► Main properties derived from this experiment: 

• Liquid drop out curve 

• Volumetric factor of the wellstream during depletion 

• Wellstream composition during depletion 


► Remarks: 

• Be sure to have vapor and liquid at each step 

• Problems imply negative calculated liquid compositions. 

• Plot Ki(P) = Yi(P) / Xi(P) 

• All lines monotonic and smooth, no crossing, ordered with volatilities 

• Use the Hoffman-Crump-Hocott plot: 

• Log(KiP) is linear regarding B(l/Tbi-1/T) with B=(log(Pci) -log(Pref)) / (1/Tbi-l/Tci) 


Fluid Studies - PVT 


!FP Training 


80 







Constant volume gas depletion 


G 


1 2 
) ■» 


(T3D 

uon 



r 


3 t 


PI =(P2at Pl^3«peat P 3 


PB 
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Constant volume depletion 




Step by step description 

0: Measure reference Volume Vsat a Pl=Psat 

1: Decrease P to P2<P1 lowering the piston, 
condensation occurs 


2: 

3: 

1: 


Reset the cell volume to Vsat, pushing up the piston 
Measure Oil and gas volumes 
Decrease P to P3<P2, ... 



| Gas 

D ° N 


Current reservoir pressure 

> 

Dew Point Pressure 

Well-stream composition 
remains unchanged 


Liquid deposit is assumed to be not mobile 


PI = Psat P2 < Psat P2 P2 P3 < P2 
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Gas unit performances 



► Range of work: 

• Pressure: 1 bar - 1500 bar 

• Temperature: -20°C to 200°C 

• Viscosity: 0.05 to 1 cPo 

► Accuracy: 

• Volumes: 

• resolution 5 10' 5 cm 3 

• accuracy 10 3 cm 3 

• Temperature: ± 1°C 

• Pressure: 

• HP: ±1 bar 

• BP: ±10' 3 bar 

• Viscosity: ±10 2 cPo 
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Liquid — Vapor equilibrium 



► Data: 


Z. =X. - L + Y. -V 

i i i 

L + V = 1 
K. = Y./X. 

i i i 

IX =ZY =1 


Balance by component 
Global balance 
Equilibrium constant 
Balance by phase 


► Calculations: 


7.. = x. •(!- v 


so : X. =Z. / 

i i 


.+ K : -X : V 

1 


Y. =K. -Z. / 1 + V 

I II 


v(k.-i)] 

i^-4 


► Rachford-Rice equation: 

f(V) = ZzjK -l|[l + v(K -l|=0 


resolution by Newton's method 

v (n+1) = v (n) — [f(v (n) )/f'(v (n) )J 
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Liquid — Vapor equilibrium 



Flash 


y Z -( K - ~ =o 

h i + v(k, - i) 


Bubble point 

£k.z,=i 


1=1 


Dew point 

n Z 
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Liquid — Vapor equilibrium 



Cal calculation of K values 



Yl 


► Experiments: 

• Binary or ternary diagrams, 

• Detailed analysis. 

► Iterative calculations: 

• Estimated initial values. 
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Liquid — Vapor equilibrium 




Yl=0.95 

Y2=0.05 

K2=0.10 
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Liquid — Vapor equilibrium 
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N2 + C 1 
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Liquid — Vapor equilibrium 



Ternary diagram 


Cl 


C5+ 
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Liquid — Vapor equilibrium 




Ternary diagram: 
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Liquid — Vapor equilibrium 

Initial value for Ki 


WILSON'S correlation 


Co = ~ln(l0) 




k i =exp C 0 •(l + w i )- 


T c 

T 


VI 


J) 


Ci 

p 
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WHITSON's correlation 


„ 69 • M- 4200 _ „ 

C i= C 2 =1 ~ 


, xO.6 

P-1 ' 


14.696 


vC,- 1 / 


f D V ’" 1 


k = 




•expjc 0 -C 2 -(l + w.)- 




M 


1 __^ 


J) 


>A 
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Liquid — Vapor equilibrium 



► State functions 

► Free enthalpy G (Gibbs) 

• G = H - T.S = (U + PV) - TS 

• dG = - S.dT + V.dP 

► At the equilibrium some of these functions display a minimum. 
In particular: dG = 0 (T, P system) 



g, = g, + RTIn 
f : fugacity 


ft 

if 


\ 

) 


Liquid - Vapourequilibrium 


f v =f L for all "i" 
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Liquid — Vapor equilibrium 



► State variables 

• Temperature, Pressure, Volume, ... 

► State functions 

• Internal energy U 

• dU = dQ + dW 

• For a reversible transformation: 

• dQ = T.dS (S: entropy) 

• Enthalpy H 

• H = U + P.V 

• dH= T.dS + V.dP 

• Free energy A (Helmholtz) 

• A = U - T.S 

• dA = - S.dT - P.dV 

• Free enthalpy G (Gibbs) 

• G = H - T.S 

• dG = - S.dT + V.dP 
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Liquid — Vapor equilibrium 



► Thermodynamics potentials 

• At equilibrium, some of these state functions display a minimum; by analogy with 
potential energy, these functions are called 

THERMODYNAMICS POTENTIALS 


► In particular: 

• dG = 0 (T, P system) 

• dA = 0 (T, V system) 
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Liquid — Vapor equilibrium 




Disturbance and return to equilibrium such that: 

dG = dG v + dG L = 0 

i.e. Y,Si dn i =~YjSi' dn f with dnf = —dnf 

i i 

Zfef -gfK =° v/ 
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Chemical potential gi and fugacity fi 



g 


i = g° + RT In 

Liquid - Vapor Equilibrium 



\ 

J 


fi v = fi L for all "i" 

Reference: Ideal gas at T and P 


fi° = XiP 

Coefficient of fugacity 

<p , = f' 


Xi P 

At equilibrium 


K i = 


y< 

Xi 


(pt 
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Equations of state 


Fluid Studies - PVT 


MARIOTTE 

1650 

p= rt 

V 

VAN DER WAALS 

1873 

0 RT 

a c 

n 

i 

> 

i 

L 

V 2 

CLAUSIUS 

1880 

D RT 

a c 

n 

i 

> 

i 

L 

T(V + c) 2 

BERTHELOT 

D — RT 

a C 

n 

i 

> 

i 

TV 2 

REDLICH KWONG 

o R T 

a c T 0 5 

n 

i 

> 

i 

v(v + b) 

WILSON 

P - RT 

A(T) 

1964 

V-b 

V(V + b) 

SOAVE REDLICH KWONG 

1972 

D _ RT 

B(T) 

n 

i 

> 

i 

L 

V(V + b) 


0 RT 

C(T) 

PENG - ROBINSON 

1976 

n 

i 

> 

ii 

L 

'n 

i 

+ 

'n 

+ 

i 

PATEL TEJA 

1981 

0 RT 

D(T) 

n 

i 

> 

ii 

L 

V(V + b) + c(V - b) 
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Acentric factor w 



1/Reduced temperature (1/Tr) 


1 1.25 1/0.7 1.5 
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Fugacity 



General expression for a pure component: 



ff^i 1 , 

(•inf C R-n 


In 

— =— — 

[ P --TT- 

dV + Z-ln Z-l 


V P J RT * 

( V v V ) 



Peng-Robinson's Equation: 

Fugacity coefficient: — 

p 

Infij = Z lLog(Z B) 

A (Z. + 2.414 -B^ 

_ 2a/2B VZ-0.414 bJ 

IFP Training I ios 
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Combination rules 



Pure components ^ mixture 

RT a m (T) 


V-b m V(V+bJ+b m (V-bJ 

b.=i>,'b, 


i=l 


a m =ZZ X ;-Xj-(l-k ij )-7v^ 


j=l i=l 


k = Binary interaction coefficients 
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Liquid — Vapor equilibrium 



Binary diagram: C02-iC4 
(experiment and EOS calculations) 



X(i)Y(i) 


Fluid Studies - PVT 


IFP Training 


110 


Fugacity 



Fugacity coefficient of component k in a mixture 


Peng-Robinson's Equation: 


r * \ 


Log 


v Px .<y 


= XT-(Z-l)-Log(Z-B) 
b_ 


r " a 

2 2>i a i( < b 


2V2B 

x Log 


a m b 

m m 

v y 

Z + 2.414B^i 


Z - 0.414B, 
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with: a |k = (l - k |k )• ^/a, • a k 
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Liquid — vapor equilibrium 
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Liquid — vapor equilibrium 



Generalized NEWTON's method 

► Target function F = f L — f v (ncomponent vector) 

• n variables 

• X 1; X 2 ,..., X n _ t and L -> R (vector of iteration variables) 

► Newton's method: j a (R a+1 -r") = f a 

(J :F(R)jacobian) 
j a =dPj8L i € [l,nc] 
j ik =ap/ax k i e [l,nc] 
k € [l,nc - 1] 

► For the next iteration, the solution vector R is given by: 

Rl + 1 =R ‘ _[/]-* 

(quadratic convergence) 
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Standard compositional grading with depth 



Variation of composition with depth 

► As hydrocarbons move into the reservoir rocks from source, they are 
influenced by the Earth's gravity, the temperature field, and capillary forces. 

► Gravity forces tend to segregate the heavier and lighter components. 

► Thus, lighter components move upwards while heavier components move 
downwards. 

► Temperature acts in an opposite direction and reduces the intensity of the 
movements. 

► The resulting effect is the decrease in solution gas oil ratio, formation volume 
factor of oil and bubble point pressure, and an increase in the viscosity and 
gasdew pressure with depth 
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Standard compositional grading with depth 
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Gas 


GOC 


Oil 


-3400 m 
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Compositional grading with depth 
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GOC transition: evolution of the phase envelopes 




Compositional grading with depth 
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Hypercritical transition: evolution of the phase envelopes 



There is no point where the bottom hole conditions are on the phase 
envelope, the fluid is never saturated and there is no GOC 
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Compositional grading with depth 



Hypercritical transition: BH density profile 



GOC transition: BH density profile 


Density (Kg/m 3 ) 
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Compositional grading with depth 




Hypercritical transition: Total GLR profile 
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Compositional grading with depth 



► Due to the gravity (y) effect, the chemical potential "g" varies with depth: 

d gj = nrij y dh 

► We know that: g ; = gj° + RT Ln (f ; / fj°) 

d gj = RT d (Ln fj) 

(nrij y / RT) dh = d (Ln fj) 
so: mj yh / RT - 0 = Ln fj h - Ln f ; 

fj h / fj = exp (mj yh RT) 
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Compositional grading with depth 

Algorithm 
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Matching Methodology 



Methodology 


► As few EOS parameters as possible 

► Most poorly defined components (highest uncertainty) 

► Maintain properties monoticity 

► Change properties by a few percent, if possible 

► Weight important measurements 

► Require "match" to reservoir and surface conditions 

► Swelling test for injection schemes? 

Trial and error until acquire "engineering feel" 
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Matching Methodology 



Prior to matching 


► Check measured data for consistency and quality 

• Compositions sum to 100%? 

• Pressure-dependent data: correct trends? 

• Material balance on CVD? 

• Property definitions? 

• Consistent units? 

• Plus fraction description? 

► EOS: Use three-parameter model - extra degree of freedom in c ; (Volume 
Shift) 
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Matching methodology 



► Choice of the variables 

• T c , p c , co of Cn+ fraction(s): saturation pressure, liquid dropout, etc. 

• Volume shift: Z-factors, densities, etc. 

• Z c or V c for LBC viscosity 

• Parachor for interfacial tension 

► Dangers 

• The "traditional" matching procedure is to adjust k cl Cn+ to obtain measured p sat 

• May imply a phase diagram deviation at low temperature. 

• Engineer must decide if deviation will affect the simulation results or separator flashes 
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Matching methodology 


Dangers 




Matching algorithm 
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Influence of Tc on the liquid deposit 
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Influence of Pc on the liquid deposit 




Pressure (bar) 
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Influence of u> on the liquid deposit 




Pressure (bar) 
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Influence of Kij on the liquid deposit 



> 

CL 

Q. 

o 


o 

Q. 

CD 

"O 

[O 

C" 



0 100 200 300 

Pressure (bar) 
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The "split" of the heavy cut 




Pressure (bar) 
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Influence of Tc on the envelope 




Temperature (°C) 
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Influence of Pc on the envelope 




Temperature (°C) 
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Influence of u> on the envelope 
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CD 

-£ 

QJ 

l . 

3 


to 

to 


CU 

l - 

Q. 


C0c cn+ + 20% 



Temperature (°C) 
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Influence of Kij on the envelope 
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ru 

-Q 


CD 

3 

Mi 

Mi 


Q) 

L- 

O. 



Temperature (°C) 
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The "split" of the heavy cut 





CD 

CD 

3 

V) 

V) 

CD 

i- 

O. 



Temperature (°C) 
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Volume shift 




v = v 


EOS 


Matching of both bottom hole and surface densities 
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Results of a 14-component match 




Measured 

Simulated 

Relative error 


values 

values 

% 

p 

r rosee 

353.000 

353.000 

0.00 

Pv a P sat 

0.315 

0.315 

0.00 

Pv a Pfond 

0.312 

0.313 

0.32 

GOR sep 

1631.800 

1733.542 

6.24 

Pl sep 

0.722 

0.730 

1.11 

GOR total 

2055.770 

2071.080 

0.74 

p L Std 

0.779 

0.767 

1.54 

Simulation of liquid deposit 

Pressure 

Measured 

Simulated 


bar abs 

values 

values 


353.0 

0.000 

0.000 


350.0 

0.135 

0.230 

See figure 

334.5 

1.111 

1.396 


318.0 

2.368 

2.596 


281.0 

5.075 

5.128 


245.5 

7.359 

7.314 


201.0 

9.574 

9.558 


159.0 

10.824 

10.952 


120.0 

11.360 

11.500 


82.0 

11.313 

11.303 


42.0 

10.722 

10.066 
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The "split" of the heavy cut 



CO 
_Q 


cu 

i_ 

3 

l f) 
IS) 

CD 



Temperature (°C) 
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The "split" of the heavy cut 



Splitting 

► "Insufficient description of heavier hydrocarbons reduces the accuracy of PVT 
predictions" Whitson C.H., SPEJ, p. 683, Aug. 1983 

► Condensates and Volatile Oils are particularly sensitive to the composition 
and properties of the plus fraction 

► Laboratories tend to give very limited analysis to the plus fraction, i.e., M N+ , 

Yn+ 

► Requires an automatic procedure to characterize the plus fraction 

► Whitson (1980) 

► Pedersen et al (1982) 

► Katz (1983) 

► Pedersen et al (1985) 

► Whitson et al (1986 & 1989) 


Fluid Studies - PVT 


!PP Training 


The "split" of the heavy cut 



Analytical Splitting 


► In splitting the heavy ends using analytical techniques, important 
requirements must be fulfilled: 

N+ N+ 

Xci = c; X ciMw i = c ^ Mw 7 

i=n i=n 

► The sum of product off the mole fractions and molecular weight divided by 
the specific gravity of each individual component is equal to that of Cn+ 


i 


Ci Mw, 
Yi 


Z 7 + MWy 

Yi 


Ahmed (1989) 

IFP Training I i48 
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The "split" of the heavy cut 



Whitson Analytical Splitting 


► Knowing the Molecular Weight and density (specify gravity) of the Con- 
traction, Whitson splitting calculation uses a three-parameter probability 
density function. 

► This procedure splits Cn+ fraction into many (30-40) small pseudo- 
components, then lumps them to the desired number. 
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4. Characterization 
and modeling 
of heavy fractions 





Heavy cut properties 



► There are so many components in petroleum reservoirs that it is impossible 
to obtain their exact composition. 

► Only pure components such as N 2 , C0 2 , H 2 S, C^-Cg up to C 45 may be identified, 
and the rest is lumped and called heavy ends. 

► These heavy ends (or plus fractions) may comprise very large number of 
isomers, and can be defined by boiling point fractions only. 

► These heavy constituents have considerable impact on the performance of 
EOS. 
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Heavy cut properties 



► For use in EOS, the components should at least be characterized by their 

• critical pressure 

• critical temperature, 

• critical volume, 

• acentric factor, 

• molecular weight, and 

• density at standard conditions. 
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Heavy cut properties 



CAVETT's method 


► Required data 

• Boiling point Teb 

• API = 141.5 /SG- 131.5 


► Tc and Teb in Fahrenheit 

T c =F 0+ FJ eb+ FJl + F 3 T eb API 


+ FJ e ;+FJ e 2 b API + F G T: b APP 



F 

G 

0 

768.071 

2.829 

1 

1.7134 

n <mi oa 1 o ' 3 

2 

- 0.10834 10 2 

- 0.30475 10 s 

3 

- 0.89213 10’ 2 

- 0.20876 10 ^ 

4 

0.38890 10' 6 

0.15141 lO -8 

5 

0.53095 10' 6 

0.11048 10' 7 

S 

0.32712 10 7 

- 0.48270 10‘ 7 

7 


0.13950 10' 9 


► Pc in psia is given by: 

Log(Pc) = G„ + G 1 T eb + GJi + G s T, b API + GJ, 3 t 
+ G 5 T e 2 b API + G 6 T eb API J + G 7 T, b API 2 
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Heavy cut properties 



Lee-Kessler 


/? c (psi) = 8.3634- 


0.0566 

Y 


0.24244+ 2.2989/ 4 


0.11857 


A 


+ 


v 

f 

1 

V 

f 

( 

V 


r 


xio ~ 3 r 


J 


1.4685 + 


3.648^1 


0.47227 


7 t 2 


xWT u 


0.42019- 


r J 

1.6977^ 

r 2 , 


r 


10 rp3 


x\0~ LU T 
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T (R) = 341.7 + 811.0Y + (0.4244 + 0.1174v)T b 
(0.4669 -3.26223y)x10 5 
+ \ 

with y - SPGR 
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Heavy cut properties 



Lee-Kessler 


f 


-In 


f Pc ^ 6.09648 

-5.92714 + + 1.28862lnT -0.169347T 


A 


co = 


14.696 


I 


b,r 


b.r 


V 


I 


b.r 


y 


(T b , < 0.8) 


co = -7.904 + 0. 1 352M» - 0.007465M; + 8.359T b r 

w d,i 


+ 


(1.408 - 0.01065M W ) 


b,r 


(T b ,r > 0.8) 


b,r 


15.2518- 15-6875 -13.4721lnT h +0.43577T 6 

b,r b,r 
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Heavy cut properties 



Standing (Mathews, Roland, Katz Correlation) 


Pc = 1188.0 - 431.0 . log (MW- 61.1) 

+ (2319.0 - 852.0 . logMW - 53.7) . (p - 0.8) 

Tc = 608.0 + 364.0 . dlog 10(MW - 71.2) 

+ (2450.0 . log(MW) - 3800.0) . dlog 10(p) 

BPR = BP + 459.67 
w = 3.0 . (log((Pc) - 2.0057) 

/(7.0 . T/BPR - 1.0)) - 1.0 
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Heavy cut properties 


RIAZI and DAU BERTs method 


► Required data 

• Boiling point (Teb) 

• Specific gravity 

► Teb in Rankine 


e= 


atd' 


(JL) = 


f \ 

iQgipPc 

T '-i 


V Teb 


/ 


0 

a 

b 

C 

M 

4.567300 10 5 

2,19620 

-1,0164 

Tc (*R) 

0.242787 10 2 

0,58848 

0,3596 

Pc (psia) 

3.128100 10’ 9 

-2,31250 

2,3201 

Vc (tf/lb) 

7.521400 10- 3 

0,28960 

-7,6666 
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Heavy cut properties 


Comparing the methods: Pc 
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Heavy cut properties 



Comparing the methods : Tc 
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Heavy cut properties 



► JOBACK's method 

P c =(0.113+0.0032- n A -£AP)" 2 

v c = 17 - 5 +Z AV 


► N A = number of atoms in the 
molecule 


Groups 

AP 

AV 

out of a cycle 



-CH3 

-0,0012 

65 

-CH2- 

0 

56 

>CH- 

0,002 

41 

>C< 

0,0043 

27 

=CH2 

-0,0028 

56 

=CH- 

-0,0006 

46 

>C= 

0,0011 

38 

=C= 

0,0028 

36 

in a cycle 



-CH2- 

0,0025 

48 

>CH- 

0,0004 

38 

>C< 

0,0061 

27 

=CH- 

0,0011 

41 

=C< 

0,0008 

32 
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Heavy cut properties 


FEDOR's method 

T c = 535 • log£ AT 


Groups 

AT 

-CH3 

1,79 

-CH2- 

1,34 

>CH- 

0,45 

>C< 

-0,22 

=CH2 

1,59 

=CH- 

1/4 

>C= 

0,89 





Viscosity 



Lohrenz-Bray-Clark's correlation 


[I(z,T ci )] 1/6 

^■[X(Z l Mw i )f 2 .e(z,P ci )r 

i i 

* (0.1023+8,0 + a,p 2 + a,p 3 + a.p 4 ) 4 

5 

a 1 =0.023364 a 2 =0.058533 
a 3 =- 0.040758 a 4 = 0.0093324 


Matching 

• Liquid: Heavy cut Vc 

• Vapor: Methane Vc 

• If necessary: a v a 2 , a 3 , a 4 


Fluid Studies - PVT 


!FP Training 


167 


Viscosity matching 




Pression (bar abs.) 
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Sampling 



Fluid Studies - PVT 



Separator 
(Tsep, Psep) 


; / 

► Gas @ Tsep, Psep 


Open hole or cased hole 
(monophasic flow 
as long as P > Psat) 


Reservoir fluid 


Water 


/rr / rcmling 


State of fluid during the production 
of a gas or oil reservoir 


172 


Sampling 


Two different groups of technique 

► Bottom hole sampling 

• Without phase separation at the sampling point: 

• open hole: RFT, MDT, RCI, ... 

• cased hole: Wireline, DST 


► Surface sampling 

• Without phase separation at the sampling point: well head, isokinetic, 

• With phase separation: at separator conditions 

RFT= Repeat Formation tester 

MDT = Modular Dynamic Formation tester 

DST = Drill Stem Test 
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Drill stem test 








- 

M 

f 

gg 

-PACKER 

* 4 ***#'■ • 

*, / * , r/T ■ . 




* : * ■ , • ■ 

i* t * * • , 

L% * * 4 *„ - . t * 




la 

b PACKER 


Flfar* 18-2 A drill stem 
test. 
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Repeat formation tester 



Pressure 

Gauge 





Pretest 

Chambers 


Equalizing 
Valve 


Sample 

Chambers 




RFT Pretest 
and Sampling Principle 

(Courtesy of Schfumberger) 
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Sampling 


► There is always mud filtrate 
invasion near the wellbore. 


► The main problem of Open Hole 
sampling is to reduce the amount of 
filtrate in the chambers 
-> Pump out module 
-> Analyzer module 



ZOne Sample 
module 


Pump-out 

module 
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The MDT tool 




Pump-out module 

Multi-sample module(s) 
(six 450cc samples) 

Sample chamber(s) 

(1, 2 3/4, or 6 gallon) 

Optical fluid 
analyser module 

Single probe 
module 


Packer module 
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Sampling 



Pressure Drawdown versus sample quality 



► Flowing well has pressure drop of Dp in near-wellbore region 

► Sample taken at points 1 & 4 are single phase 
¥ Sample taken at 2 & 3 are two-phase 
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Sampling 


Open Hole sampling 


Advantages 

► Precise depth sampling 

► Characterization of the whole fluid 
column in one shot 

► Small AP: monophasic sampling 

► Avoid doing a DST 


Drawbacks 

► Contamination by mud filtrate 

► Not representative for trace or 
aggressive components (Hg, H 2 S, ...) 

► Small sample volume, not sufficient 
for crude assay 

► Risk of emulsion for heavy crude 
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Sampling 



Surface sampling 


Reasons 

1. The difficulties in obtaining valid 
downhole samples in gas- 
condensate wells has made 
surface sampling at the separator 
the most usual method of PVT 
sample taking. 

2. If it is known in advance what the 
saturation pressure of the fluid is, 
and it happens to be below 
wellhead flowing pressure (at 
wellhead temperature) 


Fluid Studies - PVT 


Drawbacks 

► The PVT laboratory has to re-create 
the original fluid from a liquid 
sample and a gas sample!!! 

► The risks of contamination and 
leakage are doubled because of 
this. 

► The successful recombination also 
relies on accurate and repeatable 
measurements of both gas and oil 
rates during the test. 

► These measurements may require a 
revision in the light of the 
compositional analysis. 
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Sampling 


Bottom hole or surface sampling 
(main points) 


Surface 

► Metering 

► Stability (T, P, Qo, Qg, ...) 

► Separation efficiency 

► Total lift of the liquid phase 


Bottom hole 

► Well conditioning (clean up, choke, 
flowing or shut in well, ...) 

► Segregation (re opening of the well ...) 

► Small sample volume 
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Sampling 


Surface sampling: oil 




► Constant GOR as long as: 

• flowing pressure > bubble point pressure 

► Increasing GOR as soon as: 

• flowing pressure = bubble point pressure 
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Sampling 


Surface sampling: saturated oil 



► Whatever the choke, the GOR increases 

► Surface sampling is not reliable 
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Sampling 




► With the 1/8" choke, liquid accumulates at the bottom of the well: increasing GOR 

► With bigger chokes, the speed of the gas raises, the accumulated liquid is carried up 
by the gas: the GOR decreases 

► The same sequence happens again many times: zigzag curve sluggy flow 

► The flow rate must be increased (speed > 3 m/s) 
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Sampling 


Separator efficiency 


► If the separator is not 
correctly sized regarding the 
rates being put through it: 

• Carry over for gas condensate 

• Carry under for heavy / foamy 
oil 


FLOPETROL (42x10) SEPARATOR 
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Sampling 



SAMPLING 



Condensate gas compromise 

• Reduce the pressure drop 

• Lift the liquid 


Diphasic Monophasic 



Pressure 


Static 


Dew point 


Flowing 


Distance 
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Sampling 



SAMPLING 


Condensate gas compromise 

• Reduce the pressure drop 

• Lift the liquid 

A 


Pressure 



. Static 

Flowing 

Dew point 


=> Coiled tubing method 
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Sampling 


Grid to homogeneize 

velocity profile 



black / volatil 




O □ A 
□ A 
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Sampling 


n-paraffins distribution 




Comparison between MDT and representative BH sample 


it 10 


> 1 


♦ Reference sample 
Regular distribution 
■ Poluted MDT 


Mud oil ^ 


A \ 




1 * . 


20 25 30 

Carbon number 
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Analysis 



Decontamination: composition 


Decontamination 



Composition of OBM from two samples at 
the same depth 


Composition of OBM from one sample based 
on regular HC distribution 


Decontamination 
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Analysis 


Decontamination: properties 


► Volume additivity: 

• Mw 

• BH and Std fluid density 

• GOR 


• Bo 


► Extrapolation to zero pollution: 


Multi contamination 
rates 


Fluid Studies - PVT 



191 


Sampling 




Multi contamination rates 
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MDT after an OBM drilling 



Sampling 


Recommendations 



Connectivity? 
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Compositional grading? 



Sampling 



Assessment of PVT Reports 


► Sample composition check 

• Mole percentages sum to 100? 

• Nitrogen content - possible air contamination? 

• H2S concentration 

• was any value reported at the wellsite 

• are the lab measured values comparable 


► Mass balance check 

• Check reservoir fluid density and composition using process, CVD or differential data 
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Analysis 


Gas chromatography analysis 

Syringe (liquid injection) 




^r/tiS^mgEUa, 14 : 56:39 
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Analysis 




X X 


I 




_____ JMUfUr 
lit von 1 nipbttoei la c 
2 - 2-5 uviethflkmae 
r«u-sTifD< 


j.__ T 

-2-5 tri 

m-iji* 

J*J cjaetfejUtpcifit 
J-t d iiethjlbecUflt 

1 £] tl criietiyLcTclohtxtBt 
4 ictfeyl ecus* 

2 ittirl ectiDi 

j eEhjlhtptiat 
Mtlfl otciae 
Crtho-ijjlent 

1 - 1-2 uimbylcrjlobtMot 
liopraajlejelchfitnt 
1 c 3 el ctiittif 1 eye lch eiine 
1 ictbfl tl ttbflcTclohAiiDi 
tcbjl c 3 tCbfleftlcbtxaDt 


1 itchyl 

land nsfline 
ClIlH 

ICfCl 

. sbtuat 

fetbylejclcfctzint 
octane 
9 (tUC 


\ S&l il IH foie 

war*' 

! iopT]bmtat 
-E liittfcfl ocUnt 
2-3 dUttbjl acuse 
3-1 ^aeUi) setAM 
4-5 inttlfl oeciEt 
calacnt 
its* 

It A A 
anoint 

_ fl aoaiat 

1 - 3-5 tmtclylbeBitsr 
3 mbjl j ocaiot 

1 * 2-4 cr tit tbflfcenitDt 
See ba iTlbtaitnt 
ItebfltflbtAEtBC 

Ttrb«tm«itnt 

lartil iftaat 

1 - 2-3 tziiitbrlbtajtaf 


3 *Ebjl Cali 
j tcbjl call 
1*3 «tiH colt 
4 atlbjl 1 
3 itUr! 


Gas chromatography 


Fluid Studies - PVT 


!FP Training 


198 




Analysis 


Heavy fraction characterization 
GPC analysis 




V 


C$ 0 / 7 l 

c*i/m 


DO /114 

so/m 

ao/in 


ct l:i 

f r*;t jcn 
tuccjsfl 
{ lie t jcb 
lrictjoi 
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faction 
It ittiflt 
fraccian 

triJtjsa 
f riCtJ i>n 
flKtlll 


i. in 

i.tai 

J.dl 

3.711 

i:ffi 

1| 

$ 
.007 
.00 3 
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Measurements: analysis 



"Heavy" fraction 
Criteria -> Process 

Boiling point -> 

Distillation or gas chromatography 

Molecular weight -> 

Gel permeation chromatography 

Chemical family 
Liquid Chromatography 

Functional group -> 
Absorption spectrometry 
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Fluid synthesis 


Tool box 



► Samples expertise 

► Data standardization 

► Statistical analysis 

► EOS fluid model 

► Static gravitational model 

► GWD (gas show) 
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Light ends lumping methods 



C6cut: lumping according to boiling temperature (36.5<Teb<69.2 °C) 


Composants 

M 

Tc 

Pc 

w 

Teb 

% mol. 

Cyclopentane 

70,14 

511,6 

45,09 

0,1923 

322,05 

0,0887 

2-2 dimethyl butane 

86,18 

488,7 

30,80 

0,2310 

322,85 

0,0088 

2-3 dimethyl butane 

86,18 

499,9 

31,31 

0,2473 

331,15 

0,0299 

2 methyl pentane 

86,18 

497,5 

30,10 

0,2791 

333,35 

0,2629 

3 methyl pentane 

86,18 

504,5 

31,21 

0,2750 

336,35 

0,1576 

Normal hexane 

86,18 

504,4 

29,69 

0,2957 

341,25 

0,5191 

C6 (Teb) 

84,85 

503,1 

30,87 

0,2781 


1,0670 


C6cut : lumping according to carbon number 


Composants 

M 

Tc 

Pc 

w 

Teb 

% mol. 

2-2 dimethyl butane 

86,18 

488,7 

30,80 

0,2310 

322,85 

0,0088 

2-3 dimethyl butane 

86,18 

499,9 

31,31 

0,2473 

331,15 

0,0299 

2 methyl pentane 

86,18 

497,5 

30,09 

0,2791 

333,35 

0,2629 

3 methyl pentane 

86,18 

504,5 

31,21 

0,2750 

336,35 

0,1576 

Normal hexane 

86,18 

504,4 

29,69 

0,2957 

341,25 

0,5191 

Methylcyclopentane 

86,18 

532,7 

37,90 

0,2395 

344,95 

0,2628 

Benzene 

78,11 

562,1 

48,94 

0,2100 

353,25 

0,2969 

Cyclohexane 

84,16 

553,4 

40,73 

0,2144 

353,85 

0,4549 

C6(CN) 

84,25 

523,3 

35,15 

0,2521 


1,9929 


M: molecular weight 
Tc: critical temperature 
Pc: critical pressure 
CD: acentric factor 
Teb: boiling point temperature 


g/mol 

Kelvin 

bar 

dimensionless 

Kelvin 
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Calculation of mixture properties 




n n 7 Z T T 



i = 1 7=1 


% p c, 





U n U m 



z-z M M k 

i j i j ij 


i~L n j-L m 


Un U m 


MMYz.Yz. 

Tl TYl l J 

i = L n j=L m 
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Pure components and light cuts 



Pure compounds 


Components 

M 

Tc 

Pc co 

% mol. 

N2 

28.01 

126.2 

33.94 0.0400 

0.326 

C02 

44.01 

304.2 

73.76 0.2250 

2.590 

Cl 

16.04 

190.6 

46.00 0.0115 

64.685 

C2 

30.07 

305.4 

48.84 0.0908 

8.998 

C3 

44.10 

369.8 

42.46 0.1454 

5.118 

iC4 

58.12 

408.1 

36.48 0.1760 

0.737 

nC4 

58.12 

425.2 

38.00 0.1928 

2.035 


Lumping according to the carbon number 


Components 

M 

Tc 

Pc 

CO 

% mol. 

iC5 

71.91 

465.42 

34.39 

0.2231 

0.754 

nC5 

72.15 

469.60 

33.74 

0.2273 

0.902 

C6 

84.25 

623.33 

35.15 

0.2521 

1.993 

C7 

97.46 

558.34 

32.98 

0.2801 

2.048 

C8 

111.24 

586.58 

29.01 

0.3180 

1.467 

C9 

125.80 

610.23 

26.29 

0.3710 

0.993 

CIO | 

139.16 

635.02 

24.82 

0.4101 

0.808 | 


Lumping according to boiling point temperature 


Components 

M 

Tc 

Pc 

CO 

% mol. 

iC5 

72.15 

460.36 

33.84 

0.2272 

0.666 

nC5 

72.15 

469.6 

33.74 

0.2273 

0.902 

C6 

84.85 

503.13 

30.87 

0.2781 

1.067 

C7 

90.53 

538.64 

34.31 

0.2683 

1.905 

C8 

103.57 

572.13 

31.33 

0.2817 

2.084 

C9 

118.09 

601.84 

28.05 

0.3523 

1.130 

CIO 

133.72 

633.17 

26.33 

0.3851 

1.211 
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Data standardization 




GECO 

ELF 

CORELAB 

N2 

0.37 


0.35 


0.36 


C02 

2.65 


2.59 


2.67 


Cl 

66.18 


64.69 


65.69 


C2 

9.02 


9.00 


9.03 


C3 

4.93 


5.12 


4.89 


iC4 

0.72 


0.74 


0.70 


nC4 

2.01 


2.04 


1.86 


iC5 

0.67 


0.75 


0.6 


nC5 

0.95 


0.90 


0.86 


C6 

1.15 


1.99 


1.10 


C7 

1.86 


2.05 


1.84 


C8 

1.95 


1.47 


2.08 


C9 

1.15 


0.99 


1.25 


C10+ 

6.38 


7.35 


7.07 


MC10+ 

272 


236 


255 
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Measurements: analysis 



"Heavy" fraction 
Criteria -> Process 

Boiling point 

Distillation or gas chromatography 

Molecular weight 
Gel permeation chromatography 

Chemical family 
Liquid Chromatography 

Functional group -> 
Absorption spectrometry 
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Data standardization 




GECO 



ELF 



CORELAB 



N2 

0.37 


0.37 

0.35 


0.35 

0.36 


0.36 

C02 

! 2.65 


2.62 

2.59 


2.59 

2.67 

7 

2.65 

Cl 

66.18 


65.42 

64.69 


64.69 

65.69 

F 

65.06 

C2 

9.02 


8.92 

9.00 


9.00 

9.03 

F 

8.94 

C3 

4.93 


4.87 

5.12 


5.12 

4.89 

F 

4.84 

iC4 

0.72 


0.71 

0.74 


0.74 

0.70 

F 

0.69 

nC4 

2.01 

F 

1.99 

2.04 


2.04 

1.86 

7' 

1.84 

iC5 

0.67 

F 

0.66 

0.75 


0.75 

0.6 

. 

0.59 

nC5 

0.95 

F 

0.94 

0.90 


0.90 

0.86 

7 

0.85 

| C6 

1.15 

F 

1.13 

1.99 

1.07 

1.10 

F 

1.08 

Cl 

1.86 

r- 

1.84 

2.05 


1.91 

1.84 

r" 

1.93 

C8 

1.95 

r 

1.94 

1.47 


2.08 

2.08 

r 

2.13 

C9 

1.15 

r 

1.14 

0.99 


1.13 

1.25 

r 

1.27 

I C10+ 

6.38 

r 

7.46 

7.35 

7.76 

7.07 

r 

7.76 

MC10+ 

272 

236 

236 

236 

255 

236 


Same lumping method - same Mw for the heavy cut (C10+) 
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Liquid deposit 



o 

> 

\p 

ON 


LO 

O 
Q_ 
C U 
"O 
T3 

cr 



Pressure (bar abs.) 


1: CM E (V @ P init) 
2: CM E (V @ P dew) 


3: CVD (V @ P dew) 
4: CM E (V @ P) 
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Fluid Synthesis 



Field cases 
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7. Gas injection 


Specific PVT 
experiments 




Gas injection 




Mechanisms 

► According to gas and oil composition 
also reservoir conditions: 

• No miscibility: one phase -> swelling 
effect depends on the undersaturated oil 
in the reservoir. 



► According to gas and oil composition 
also reservoir conditions: 

► No miscibility: two phases sweeping 
efficiency depends on mobility ratio. 



Progressive change from gas to oil 



► According to gas and oil composition 
also reservoir conditions: 

• Miscibility is possible: obtaining of a 
progressive change -> means Criticality. 
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Gas injection 



► Oil: 40°API, BPP = 235 b, BHP = 360 b 


► First matching on standard oil only 
measurements 

► Experiment: gas separator injection in 
oil for several steps, bubble points 
measured 

• The model based on standard 
measurements is wrong 

• Phenomena are not well represented (L/V, 
composition, ...) 

• Necessity to take into account the effect of 
gas injection. 

► Need for gas injection experiments, 
which represent the right phenomena: 
Miscibility and Criticality 


Saturation pressure versus gas mole injected 



Gas mole injected for 1 mole of liquid 
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Gas injection 



Specific PVT experiment 
► Swelling test 

Gas is added to the initial reservoir fluid, 
mixture is kept monophasic (increasing of exp. 
pressure) 

• Saturation pressure profile 

• Criticality location 

• Relative volume (swelling effect) 


Saturation pressure versus gas mole injected 



Relative volume versus gas mole injected 


^>No phase separation, no material balance 
^Approach of criticality location 
^Standard PVT for each mixture 
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Specific PVT experiment 

► Multiple contact test (MCT) 

Gas is added at constant pressure: equilibrium 

• One phase is withdrawn (gas or liquid) 

• Compositional profile 

• Liquid relative volume (ultimate vaporization) 


Relative Liquid Volume versus number of moles injected 



► Represent the reservoir around injector 

^Experimental problems: equilibrium should be reached, phase separation V/L, 
Mat. Balance 
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Miscibility 
Slim tube 
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Miscibility conditions: MCMP 


100 JT 


90 


a so 


70 



~T 

100 

Exa-Plans 12/96 - UAA 96/S/92 


300 Pressure 
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Miscibility 




Ternary diagram 




Miscibility 



First Contact Miscibility (FCMP) 



Miscibility 



Multiple Contact Miscibility 



"Front" miscibility (vaporizing gas drive) 


The front gas vaporizes the oil intermediate components: gas becomes heavier & miscible 
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Miscibility 



Multiple Contact Miscibility 



"Rear" miscibility (condensing gas drive) 


The rear gas intermediate components condense in the oil: oil becomes lighter & miscible 
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Miscibility 



Fluid Studies - PVT 



Multiple Contact Miscibility 
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Miscible model 


► 


► 


Interfacial tension: 


The following are modeled: 



= F{cr)K rogmn+ {\-F(a))K 

= F (^_ + (i -Ha))*,. 



with F (a) = 


a V/ " 


V a o J 


(Tq « l dyne / cm 


n = 2 


to be adjusted according to the situation 
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Miscible model 




¥ The flowrates are therefore proportional to the saturations 
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Interfacial tension 


Sugden and Mac Leod's method (mixtures) 


cr = 






4 


X:Y: Molar fraction of component i in the liquid and vapor phases 


p L ,p v : Liquid and vapor phase densities 
M l ,M v : Molecular weight of liquid and vapor phases 
Pa i: Parachor of component i 
ailnterfacial tension 


g/cm 3 

g/mole 

(dyne/cm) 1/4 . cm 3 /mole 
dyne/cm 
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Matching of interfacial tension 




Pressure (bar abs.) 
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Black-Oil model 



Definitions: 



Black-Oil model 
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Black-Oil model 






Gas properties 
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Black-Oil: correlations 



Fluid Studies - PVT 


ESTIMATION OF THE BUBBLE POINT PRESSURE Pb 


to, solution GQfi {m 1 ltd / m* s*3 

15 20 30 4G 60 80 100 150 200 


Saturation pressure 



: 90°C 


Examp e: 

Temperature 
Stock tank oil density : 0.85 
Solution GOR : 1 20 m 3 std/m 3 st 

Gas density : 0.90 {sld cond.) 

2 oo Desired biible pressure : 1 65 bar 


( STANDING'S CORRELATION ) 
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Black-Oil: correlations 
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Formation volume factor estimation 


Example : 

TrouverB^j 

On donne : Rsb ■ 300 cf/Bbl 
T » 180°F 
dg-0.7 
do =36°APi 

1 Entrer Rsb * 300 (eck gauche) 

2 Vers la droite dg ® 07 

3 Descendre do = 35 

4 Vers la droite T * 1 80 

5 Descendre: Bob® 1.20 
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Differential vaporisation 


Reservoir temperature 


o 



p 2 <Pi 


T = 15 °C 



P 3 < P 2 P = 1 bar 
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Composite Bo 
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Condensate content 
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Composite black-oil quantities 



► From a PVT study: 

• process test 

• differential liberation 

1. P > P sat 

• Boc (P) = V(P) / V(Psat) * Bop (Psat) 

• Rsc (P) = Rsp (Psat) 

2. P < P sat 

• Boc (P) = Bod (P) * (Bop (Psat) / Bod (Psat)) 

• Rsc (P) = Rsp (Psat) - ARsd(P) * (Bop (Psat)/Bod (Psat)) 

• ARsd(P) = Rsd (Psat) - Rsd (P) 
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Summary 
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Compositional 

data 


Fluid Studies - PVT 


Compositional data: characterization of cuts 



COMPOSITION AND PHYSICAL PROPERTIES 

Name 

ZF % 

Mw 

Tc (K) 

Pc (bar) 

Acent fact 

Vc (cm3 ) 

Parachor 

GN1 

71.61 

16.12 

190.31 

45.88 

0.0116 

98.9 

76.6995 

GN2 

11.19 

31.28 

301.25 

49.65 

0.0957 

143.0 

105.4330 

GN3 

7.11 

48.95 

387.44 

40.49 

0.1600 

221.5 

162.1903 

GN4 

3.87 

84.05 

506.62 

32.49 

0.2621 

351.9 

268.0716 

GN5 

2.27 

122.66 

618.15 

26.34 

0.3000 

511.8 

422.4642 

GN6 

3.11 

204.81 

702.45 

18.50 

0.6000 

975.0 

481.6625 

GN7 

0.84 

340.00 

787.25 

17.50 

0.9200 

1320.0 

799.5959 
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BINARY INTERACTION COEFFICIENTS 


GN1 

GN2 

GN3 

GN4 

GN5 

GN6 

GN7 

GN1 


I -0.00155 

0.00032 

0.00890 

0.01851 

0.05840 

0.06738 

GN2 

-0.00155 


0.00840 

0.00137 

0.00278 

0.04802 

0.04800 

GN3 

0.00032 

0.00840 


0.00146 

0.00276 

0.02998 

0.02999 

GN4 

0.00890 

0.00137 

0.00146 


0 

0 

0 

GN5 

0.01851 

0.00278 

0.00276 

0 


0 

0 

GN6 

0.05840 

0.04802 

0.02998 

0 

0 


0 

GN7 

0.06738 

0.04800 

0.02999 

0 

0 

0 
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Why lumping? 


► Because of the high number of components, the analysis of a large PVT report 
may reveal impractical on a small or overloaded machine 

► The cost of EOS calculations are directly related to the number of 
components (N c * 1 2 ) used to characterize the reservoir fluid studied. 

► In reservoir simulation of the compositional type, the EOS calculations must 
be repeated numerously in each grid block and at each time step, (reducing 
the number of components -> reducing the number of equations) 

► For anything but the simplest (or economically limited) reservoir cases, a 
complete compositional simulation using an extended fluid analysis becomes 
prohibitive in time and cost terms. 

► What is required is a procedure, which will maintain the thermodynamic 
consistency while using a minimum number of representative components or 
'pseudo-components'. 
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How to lump 



According to the constraints due to 

1. Simulation: 

• Kind of fluid 

• Production scheme (gas injection ...) 

• Thermodynamic phenomena 

2. Process: 

• Surface treatment (particularly for the gas phase) 

• Specifications 
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Lumping 


Rules for grouping 

► K values with the same magnitude 

• Cj - n 2 

• C 2 - C0 2 - C 3 

• c 4 -c 5 

c - c 

• c 6 -c 10 < J? 

• C n P 8 10 


► Similarity of properties 

► Insensitivity of experiments to trial grouping 


► Example: 

• Obvious candidates 

• iC 4 and nC 4 C 4 

• iC 5 and nC 5 C 5 


Fluid Studies - PVT 


!PP Training 


257 


Lumping 



► Lumping scheme (SPE 13119) 

• Properties chosen 




m 



VO . 45724 R T 

-n/ P ci 

0 . 0778 R T , 


f( co ,) and M 
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Lumping algorithm 


Nc class centers are 
taken at random 



Each comp, is assigned to 
the closest class center 





Operation is renewed until Barycenter of each class 

no more change occurs is taken as new center 
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Properties calculation 



n n 7 Z T T 



U n U m 


r'=l JPf Pf 


c 


i=l j =x \ rczc, 

rji V J 

- Wc, 


k = 


2 - 

j — Lj n 


s 


run 


i = 1 P* 


U n u m 

MMYz.Yz. 

n n i j 

* = L n j = L m 


R = 


T, 


n 


C 


C 


” ZJc. 


2 p 

i = 1 r C t 


CO = Y J Z , 0J , 

i = 1 


n 


M = Y J X r M i 


i=\ 
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Results of a 7-component match 



Fluid Studies - PVT 



Valeurs 

mesurees 

Valeurs 

simulees 

Erreur 
relative % 

p 

~rosee 

353.000 

353.000 

0.00 

Pv & Psat 

0.315 

0.314 

0.32 

Pv ^ Pfond 

0.312 

0.312 

0.00 

GOR sep 

1631.800 

1783.837 

9.31 

Pl sep 

0.722 

0.734 

1.66 

p L Std 

0.779 

0.786 

0.90 

Simulation du depot liquide 

Pression bar 

Valeurs 

Valeurs 


abs 

mesurees 

simulees 


353.0 

0.000 

0.000 


350.0 

0.135 

0.259 

Voir figure 

334.5 

1.111 

1.553 


318.0 

2.368 

2.847 


281.0 

5.075 

5.459 


245.5 

7.359 

7.585 


201.0 

9.574 

9.654 


159.0 

10.824 

10.899 


120.0 

11.360 

11.381 


82.0 

11.313 

11.186 


42.0 

10.722 

10.008 
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Results of a 7-component match 




Pressure (bar abs.) 
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Results of a 14 and a 7-component match 




o 


Fluid Studies - PVT 


100 200 300 

Pressure (bar abs.) 


400 
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Results for different groupings 



Cut 

Mole 

Cl 

0.334 

N2 

0.007 

C02 

0.006 

C2 

0.05 

C3 

0.055 

IC4 

0.021 

C4 

0.039 

IC5 

0.02 



C5 

0.028 

C6 

0.044 

PC71 

0.096 

PC72 

0.099 



PC73 

0.087 

PC74 

0.068 

PC75 

0.046 
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U.UDD 

IC4C4 

0.06 

IC5C5 

0.048 

C6P 

0.44 



!FA 


Coupe 

IV 

flole 


C1N2 

0.341 


C02C2C3 

0.111 


IC4C4 

0.06 


IC5C5 

0.048 


C6PC71 

0.14 


PC7273 

0.186 


PC7475 

0.114 



Properties profile vs depth 



N° 

Cote en 

m/TR 

P.fond bar 

abs. 

T. fond °C 

Vis. fond 

cPo 

p fond 
Kg/m 3 

P. sat bar 

abs. 

Vis. sat cPo 

p sat Kg/m 3 

Rs m 3 /m 3 

Bo m 3 /m 3 

p std Kg/m 3 

1 

-3040.0 

390.30 

136.30 

0.0420 

320.10 

352.48 

0.0398 

302.50 

2053.38 

8.1135 

779.3 

2 

-3060.0 

390.90 

136.80 

0.0426 

323.00 

356.06 

0.0404 

306.90 

1976.57 

7.8420 

780.1 

3 

-3080.0 

391.60 

137.30 

0.0431 

325.80 

358.02 

0.0411 

310.40 

1903.34 

7.5692 

780.1 

4 

-3100.0 

392.20 

137.80 

0.0437 

328.90 

361.03 

0.0422 

314.80 

1830.09 

7.3038 

780.6 

5 

-3125.8 

393.10 

138.40 

0.0446 

333.50 

364.57 

0.0435 

320.70 

1732.10 

6.9467 

781.1 

6 

-3151.6 

393.90 

139.10 

0.0457 

338.80 

368.53 

0.0450 

327.50 

1629.21 

6.5746 

781.8 

7 

-3177.4 

394.80 

139.70 

0.0470 

345.00 

372.72 

0.0472 

335.20 

1523.92 

6.1910 

782.6 

8 

-3203.2 

395.70 

140.40 

0.0488 

353.10 

378.43 

0.0502 

345.60 

1401.78 

5.7477 

784.0 

9 

-3229.0 

396.60 

141.00 

0.0515 

364.10 

384.18 

0.0519 

358.80 

1261.60 

5.2380 

785.4 

10 

-3239.2 

396.96 

141.25 

0.0530 

370.00 

386.89 

0.0542 

365.80 

1196.22 

5.0008 

786.2 

11 

-3249.4 

397.33 

141.51 

0.0551 

377.60 

389.74 

0.0580 

374.40 

1119.60 

4.7238 

787.1 

12 

-3259.6 

397.70 

141.80 

0.0585 

386.00 

393.59 

0.0990 

387.40 

1018.49 

4.3588 

788.7 

13 

-3269.8 

398.10 

142.02 

0.0996 

468.70 

395.93 

0.1142 

467.90 

582.81 

2.8162 

799.7 

14 

-3280.0 

398.57 

142.27 

0.1165 

486.80 

391.67 

0.1246 

384.60 

519.96 

2.6020 

802.5 

15 

-3290.2 

399.06 

142.53 

0.1286 

397.30 

388.45 

0.1334 

493.90 

487.17 

2.4918 

804.2 

16 

-3300.4 

399.56 

142.78 

0.1389 

505.10 

386.00 

0.1412 

500.80 

464.42 

2.4159 

805.6 

17 

-3310.6 

400.06 

143.04 

0.1484 

511.50 

383.73 

0.1484 

506.40 

446.53 

2.3566 

806.8 

18 

-3320.8 

400.57 

143.29 

0.1574 

517.00 

381.62 

0.1553 

511.20 

431.80 

2.3081 

807.8 

19 

-3331.0 

401.09 

143.55 

0.1661 

521.90 

379.64 

0.1615 

515.40 

419.05 

2.2664 

808.7 

20 

-3341.2 

401.61 

143.80 

0.1742 

526.10 

377.81 

0.1678 

519.00 

408.38 

2.2316 

809.5 

21 

-3351.4 

402.14 

144.06 

0.1823 

530.10 

376.13 

0.1742 

522.40 

398.64 

2.2000 

810.4 

22 

-3361.6 

402.67 

144.31 

0.1907 

533.90 

374.60 

0.1799 

525.60 

389.70 

2.1710 

811.2 

23 

-3371.8 

403.20 

144.57 

0.1984 

537.20 

372.89 

0.1855 

528.40 

381.67 

2.1452 

811.8 

24 

-3382.0 

403.74 

144.82 

0.2062 

540.40 

371.40 

0.1913 

531.00 

374.56 

2.1224 

812.5 

25 

-3392.2 

404.28 

145.08 

0.2142 

543.50 

369.85 

0.1971 

533.60 

367.46 

2.0997 

813.1 

26 

-3402.4 

404.83 

145.33 

0.2223 

546.40 

368.54 

0.2026 

636.00 

361.11 

2.0793 

813.8 

27 

-3412.6 

405.37 

145.59 

0.2300 

549.00 

367.18 

0.2081 

538.20 

355.25 

2.0607 

814.4 

28 

-3422.8 

405.92 

145.84 

0.2380 

551.70 

365.79 

0.2137 

540.40 

349.51 

2.0425 

814.9 
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Pseudo process 





Sale Gas 


Fuel Gas 

► 


Condensates 



99% 


1 % 
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Pseudo process 


► Split factors are used in STAB to obtain a liquid without methane 

► The following split factors are used to restore the composition of the 
condensate over several years, 

• N2: 100%, C02: 100%, Cl: 100%, 

• C2: 96%, C3: 58%, iC4: 13%, nC4: 10% 

• iC5: 0, nC5: 0, 


Fluid Studies - PVT 


!PP Training 


267 


Lumping / ungrouping 



Hypothesis ■ are zero 

Oi L = Ln [fi L /(xi.P)] = Co + Ci (ai) 1/2 +C L 2 bi 
3>i v = Ln [fi V /(yi.P)] = Co + Ci (ai) 1/2 +d 2 b \ 

At equilibrium ■ f L =f' V 
Ln Ki = Ln (ys / x i) = O i L - O i V 
Ln Kj = A Co + A Ci (as ) 1/2 + A C 2 b i 


Ln Ki = Co + Ci (ai) 1/2 + C 2 bi 
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Fugacity 


► Fugacity coefficient of component k in a mixture (Peng-Robinson's EOS) 


Log 


f, ^ 


K p \J 


= -^(Z-l)-Log(Z-B) 
b 


n \ 

2 Z x , a ik h 


2 V 2 B 


m 


m 


V 


J 


xLog 


rz + 2.414B^ 

• i f 

VZ -0.414B J 

with: a., = / 
ik V 


a. - a. 

1 k 


when: 6., =0 
ik 
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Lumping / ungrouping 



► From: 

• 7 pseudo-components (z,) to 16 components (Z,) 

► Equilibrium @ 7 pseudo-components: 

• We know: x i; y, and 0 ^ k, 

• There are: 7 equations Ln(kj) = C 0 + C 2 (a^ 2 + C 2 bj ^ C 0 , C x and C 2 

► Hypothesis: 

• 0 is the same for the 2 representations (7 and 16 comp.). 

• For the 16 components we know: Z j; 0, C 0/ Cl and C 2 so Z t , 0 and Kj 


X, = Z, / (1 + 0 (K, -1)) 
Y| = Zj K| / (1 + © (K, -1)) 
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Lumping / ungrouping 



► For a condensate gas, a CVD can be numerically performed with a detailed 
fluid representation. 

► So, the detailed compositional profile of the produced gas versus the 
reservoir pressure is obtained. 

► At a given reservoir pressure, each pseudo component is ungrouped 
according to the new "true" proportions. 
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Central Area [phase 1] 

Eastern Area [appraised] 

Northern Area [to be appraised] 


ALKHALIJ 


Fluid synthesis 

IFPTraii ling 


Okm S 

1 \ 
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ff l 

Rostam 


Field case / A 


Properties 

Initial 

A3 DST1 

1992 

ANE 

MDT 

2000 

BH conditions 



(52.7X- 131.7b) 



Density (kg/m3) 

822,3 

833,6 

Viscosity (cPo) 

2,39 

4,27 

Saturation pressure (bar) 

65,4 

32,6 

Density (kg/m3) 

814,6 

826,9 

Viscosity (cPo) 

2,12 

3,48 

Process test 



(40°C-6. 7b / 40° C- lb / 15°C-lb) 



Total GOR 

43,2 


Boi 

1,1218 


Liquid density (kg/m3) @ 40°C 

862,1 


Stock tank liquid density (kg/m3) 

874,3 


Composition test 



(50°C-lb / 15°C-lb) 



Total GOR 


20,1 

Boi 


1,1018 

Liquid density (kg/m3) @ 50°C 


868,8 

Stock tank liquid density (kg/m3) 


891 
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Components 

Initial 

A3 DST1 

1992 

ANE 

MDT 

2000 

N2 

0,23% 

0,12% 

C02 

0,67% 

0,16% 

Cl 

17,74% 

6,91% 

C2 

4,48% 

3,49% 

C3 

4,98% 

4,65% 

IC4 

1,24% 

1,27% 

NC4 

3,27% 

3,40% 

IC5 

1,97% 

2,26% 

NC5 

2,18% 

2,48% 

C6 

4,47% 

5,35% 

C7 

4,30% 

5,57% 

C8 

4,03% 

5,05% 

C9 

3,68% 

4,97% 

CIO 

3,59% 

4,63% 

CUP 

43,17% 

49,69% 
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Field case / A 



2 wells: A3 and A-NE 


Very different compositions, 
very different PVT behaviors, 
about 10km apart with a 
possible sealing fault between them 

2 different fluids: 
compartmentalization 
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Field case / A 



Data from drilling 


Problems during drilling, 

the well has remained "under water" for 2 months 





2 different fluids: 
compartmentalization 


Washing by the mud 

filtrate 
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Field case / A 



Geochemistry study 


same origin, 
same maturity, 
no biodegradation 


i 

Washing by the mud 

filtrate 
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Thermodynamic study 


fluid model for A3, 

simulation of a multi-contact experiment 
between A3 fluid and water 


Soreide & Withson EOS 
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Field case / A 




Initial 

After water 

ANE 




Properties 

A3 DST1 

washing 

MDT 





1992 

(34 vol) 

2000 




BH conditions 








(52.7°C- 131.7b) 






After water 

A NE 






Comoonents 

washing 

MDT 

Density (kg/m3) 

822,3 

844,8 

833,6 



(34 vol) 

2000 





Viscosity (cPo) 

2,39 

3,08 

4,27 


N2 

0,05% 

0,12% 






C02 

0,00% 

0,16% 

Saturation pressure (bar) 

65,4 

24,0 

32,6 


Cl 

6,87% 

6,91% 

Density (kg/m3) 

814,6 

833,9 

826,9 


C2 

3,54% 

3,49% 

Viscosity (cPo) 

2,12 

2,60 

3,48 


C3 

5,31% 

4,65% 






IC4 

1,44% 

1,27% 

Process test 





NC4 

3,74% 

3,40% 

(40°C-6. 7b/ 40° C-lb / 15°C-lb) 





IC5 

2,29% 

2,26% 






NC5 

2,55% 

2,48% 






C6 

5,22% 

5,35% 

Total GOR 

43,2 

20,2 



C7 

5,03% 

5,57% 

Boi 

1,1218 

1,0668 



C8 

4,72% 

5,05% 

Liquid density (kg/m3) @ 40°C 

862,1 

861,2 



C9 

4,31% 

4,97% 

Stock tank liquid density (kg/m3) 

874,3 

873,5 



CIO 

4,21% 

4,63% 






CUP 

50,59% 

49,69% 

Composition test 








(50° C-lb/ 15° C-lb) 







Total GOR 


25,4 

20,1 




Boi 


1,0890 

1,1018 




Liquid density (kg/m3) @ 50°C 


861,3 

868,8 




Stock tank liquid density (kg/m3) 


878,3 

891 
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Field case / A 



From MDT measurement, a slight trend can be observed on Cl content, the 
saturation pressure and GOR as a function of the pumping duration before 
sampling 


MDT number 

Pumping 

duration 

Cl content 

(mol %) 

BH density 

(kg/m3) 

Sat. Press. 

(abs. bar) 

GOR 

(Sm3/Sm3) 

ST liq. Dens. 

(kg/m3) 

Water 

content 

90 

2h 30 

6,5 

837,2 

30,4 

19,7 

891 

50% 

128 

3h 

6,9 

833,6 

32,6 

20,1 

891 

37,50% 

129 

3h 

7,2 

833,7 

32,9 

20,1 

894,1 

40% 


No conclusion was drawn from this single trend 
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Field case / A 



Conclusions 


► The unexpected A-NE sample composition can be explained by a water 
washing process 

► A representative sample was obtained after a longer pumping period 
duration 


Later, the production data confirmed that the two fluids were identical 

(no compartmentalization) 
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Field case / V 




North fane 


as. 


Sliqhty modified fluids 


Infer mediate 
fluids 
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Field case / V 



► GOR measurements obtained from a lot of wells do not show any obvious 
trend according to perf depth and/or well location 

► Highly undersaturated, the fluid exhibits a non classical GOR profile vs. depth: 
from 280 to 150 for only 70 m depth variation! 


i' 

Different fluids: compartmentalization 
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Field case / V 



Data from Geology - Geochemistry 


2D Temis simulation on the Central Graben: 
accumulation conditions (P, T) 

(at the beginning (15 M years), saturated fluid with a gas-cap) 
Sedimentation speeding up during the last 5 M years 

Paleo gas-cap dissolution 
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Field case / V 



Data from Geology - Geochemistry 
Accumulation pressure # 250 b 

Data from Seismic 

Due to a gas "cloud" over the cap-rock 
It was difficult to interpret seismic data 

Data from PVT study 
BHP # 460 b 

BPP: top 330 b, top -70 m # 270 b 
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Paleo gas-cap dissolution 

IFP Training I 286 

EP - Fluid synthesis — ' 


Field case / V 
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Paleo gas-cap dissolution 
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Field case / V 



North 


South 



Paleo gas-cap dissolution 
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Field case / V 



Today, the calculated diffusion mean distance 
is about 500 m 


-> only a part of the reservoir has been highly affected 
by the dissolution / diffusion phenomenon 


-> original fluid (not affected by the phenomenon) 
should be produced by the flank wells 


EP - Fluid synthesis 


IFP Training 


289 


Field case / V 



Thermodynamic study 


A differential liberation of a top fluid has been performed 
from bottom pressure to 250 b (accumulation pressure) 


New calculation of the composition of the fluid obtained according 
to the depth in order to compare with the fluid flank 


Gravitational model 
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Field case / V 



Component 

Composition 

top fluid 

Diff. Lib. @ 90°C 

liquid phase (step 250 b) 

Composition 

37.5 m below 

@ 90°C 

2/8 

flank fluid 

Composition 

68 m below 

@ 91. 1°C 

2/11 

flank fluid 

N2 

0,66 

0,48 

0,48 

0,26 

0,47 

0,36 

C02 

0,35 

0,35 

0,35 

0,10 

0,35 

0,35 

Cl 

54,51 

46,40 

46,05 

44,86 

45,77 

45,68 

C2 

7,88 

7,85 

7,82 

6,29 

7,80 

6,16 

C3 

5,40 

5,76 

5,75 

5,92 

5,73 

6,17 

IC4 

0,76 

0,84 

0,84 

0,88 

0,83 

0,83 

NC4 

2,81 

3,16 

3,15 

3,54 

3,15 

3,50 

IC5 

1,00 

1,15 

1,15 

1,34 

1,15 

1,54 

NC5 

1,55 

1,80 

1,80 

1,89 

1,80 

2,17 

C6 

1,91 

2,26 

2,26 

2,31 

2,26 

2,88 

C7P 

23,17 

29,94 

30,35 

32,61 

30,68 

30,36 
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Conclusions 


Field case / V 


The classification of the GOR values 
can be explained according to: 


the distance to the anticline axis, 
the depth, 

the fracturation level 
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Field case / V 




Field case / G 



Deep offshore field 
PVT data 
MDT pressures, 

composition profile with depth 

* 

Continuous oil column 


EP - Fluid synthesis 


!PP Training 


294 


Field case / G 



Non classical evolution with depth of some parameters 
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Field case / G 



Non classical evolution with depth of some parameters 


Saturation pressure (b) 

250 255 260 265 270 275 
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Field case / G 



Data from geochemistry 

very mature underlying source rock 
producing gas 

Data from GWD 
non sealing cap-rock (Cl leakage) 

External gas flux 
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Field case / G 



a. 

LU 

O 



Gas flux from 
Source Rock 


Gas flux from 
biodegradation 
activity 


AND/OR 


Gas flux from 
intermediate reservoirs 


Source Rock 
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Field case / G 



Gas flux 
in oil pool 



Gas diffusion Oil stripped and 

Pseudo gas-cap through cap-rock saturated at the top 



Step 1 


Step 2 


Step 3 


Step 4 
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Field case / G 


Saturation pressure (b) 

250 255 260 265 270 275 

2300 


2400 


i/i 

— 2500 

-C 

+-> 

a. 

ai 

a 

2600 


2700 

The non classical evolution has been explained and 
simulated by means of a partial gravitational model 

IFP Training I 300 
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Fluid studies PVT 
(Sampling - PVT) 

Quick sum up 


/FPTraining 


Reservoir fluid classification 



Oil system : T res < Tc 

Gas system : T res > Tc 

+ Initial reservoir Tres, Pres 
★ Surface separator Ts, Ps 
■fo Critical point Tc, Pc 



Under saturated 



Temperature 



Temperature 



Fluid studies PVT 
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Standard compositional grading with depth 



Critical transition without any GOC 




3200 m 

3250 m 



Due to the displacement of the critical points, T res 
(^constant) which is higher than the corresponding Tc 
at the top, becomes lower at the bottom ^ 
continuous change from a gas phase to an oil phase 

without any interface (no GOC). 

None of the initial bottom hole conditions enters the 
phase envelop ^ the reservoir fluid is always under 
monophasic state. 
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Sampling (1) 



Cased Hole sampling 

(Sampling after casing - during Testing (DST)) 

• Surface Sampling: 

• Well Head (monophasic) 

• Specific manifold 

• Separator 

• Downhole Sampling: 

• Wireline 

• Chamber included in test assembly 


Diphasic 


Monophasic 


Open Hole sampling 

(Sampling before casing - while drilling - no flow at surface) 

• Downhole Sampling: I — S Monophasic 

• MDT type (Schlumberger), RDT (Halliburton), RCI (Baker) 



Downhole Sampling : the fluid is monophasic at sampling point, in open hole sampling the 

fluid may be polluted (after OBM drilling) 

Sampling at surface : the fluid is generally diphasic 


Fluid studies PVT 
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Sampling (2) 



Reliable sampling 


Whatever the technique: bottom hole flowing pressure > saturation pressure 


► 


Bottom Hole Sampling (BHS): 


• MDT type: fluid contamination (after an Oil Based Mud drilling) ^ suitable corrections 
to get the properties of the actual reservoir fluid 

• DST: sampling tool location w.r.t perforations (segregation, liq-vap equilibrium) 


Surface Sampling (separator): 


• Constant GOR: 

• constant separator conditions (P, T) 

• lift up velocity for wet gases (v > 2-3 m/s to ensure the total liquid lift up) 

• Heater prior the separator input (to avoid wax deposit) 

• Separator efficiency (carry-over (wet gases), carry-under (heavy oils)) 
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Sampling (3) 



Open hole sampling after an OBM drilling 




Invaded 

zone 



► Fluid sample (MDT type) is contaminated by the 
hydrocarbons added to the mud, 

► Whatever the pumping duration (for clean out), zero 
pollution is generally impossible, 

► So, PVT experiments will be carried out on a 
contaminated fluid. 


Do not use raw data from PVT study without 
any appropriate corrections 



Zero pollution 



In order to determine the whole set of 
properties of the actual reservoir fluid 


Fluid studies PVT 




At a given depth, take 3 samples at different pollution rates. 
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Saturated 







PVT experiments for oil 

► Constant mass expansion 

• bubble point pressure 



Pl>Psat P2>Psat P3>Psat P4<Psat P5<Psat P6<Psat 


► Multi-stage separation 

• initial oil composition 

• initial oil density 

• initial oil process Rs 

• stock Tank Oil (STO) properties 


- Gas treatment 



U 


TT* Fue| § as 

- Injection 


bottom 
(Initial oil) 


Reservoir temperature 


T = 15 X 


Initial 

Oil 



r 


r 


p 2 <p. 
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► Differential vaporization fluid properties vs Pressure 

• differential Bo function of P 

• differential Rs function of P 

• STO properties 
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PVT experiments for gas condensate 


► Constant mass expansion (CME) 







P > Psat PI = Psat P2 < Psat P3 P4 P5 = Psat P6 

► Constant volume depletion (CVD) 

• very close to the actual behavior 

• produced gas composition function of P 

ti 



II T 






PI = Psat P2 < Psat P2 

Fluid studies PVT 


P2 P3 < P2 P3 


Liquid deposit in the reservoir 
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Gas Injection 



Mechanisms 




Change of So ( 71 ) and jlx (^1) 


According to gas and oil composition also reservoir conditions: 

No miscibility: one phase -> swelling effect depends on the 
undersaturated oil in the reservoir 


According to gas and oil composition and also reservoir 
conditions: 

No miscibility: two phases ->the sweeping efficiency depends 
on the mobility ratio 

Gas injector: C 02 , MUM 

separator gas... 

Oil producer 




Front or fingerings 


Progressive change from gas to oil 

Fluid studies PVT 


According to gas and oil composition and also reservoir 
conditions: 

Miscibility is possible: obtaining of a progressive change -> 
means Criticality 


/PP Training 311 


Definitions (1) 


Surface 

(separators conditions) 


T i p i 


1 


Std Gas 1 


condensates 


1 bar / 15° C 


Std Gas 2 


0 

t 2 p. 


2 p 2 1 bar/ 15°C 


Fluid 


Bottom 


(Reservoir conditions) 


V oil 


Pi 


P 1 < P sat 
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When the oil and gas reservoir 
phases run through surface facilities 

For oil reservoir 

Bo = Voil/STO 
Rs = Std Gas 2 / STO 

For gas reservoir 

Bg' = V gas / Std Gas 1 

Rv = cgr = condensates / Std Gas 1 

Bg = V gas / V gas in Std cond 


► Volumes at surface (lb/15°C) depend on the process 
scheme; so, Bo, Rs, Bg' and Rv depend on the process. 

► Bg does not depend on the process. 

► Except for dry gas, Bg is theoretical. 

• Bg = (Z * T) / (288.15 * P) T in K, Pin bar 

• Bg ~ 1 / P Pin bar 

► As the reservoir pressure changes, oil and gas phases in 
the reservoir change too; so, Bo, Rs, Bg', Rv and Bg are 
functions of the reservoir pressure. 
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Definitions (2) 



Surface 

(separators conditions) 


Fluid 


i — 

■ 

t i p i 


■ 


Std Gas 1 


o - 

T 2 P 2 


condensates 


1 bar /15°C 


Std Gas 2 


0 

T, P, 


°l — ' 

T 2 p 2 : 1 bar / 15° C 


Bottom 


(Reservoir conditions) 


V oil 


Pi 


P 1 < P sat 


Gas properties 


For gas reservoir 

Bg' = V gas / Std Gas 1 
Rv = condensates / Std Gas 1 
Bg = V gas / V gas in Std cond 




For oil reservoir 

Bo = V oil / STO 
Rs = Std Gas 2 / STO 
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Calculation of STOMP 



How to calculate the Stock Tank Oil Initially In Place? 


► The geologist knows: 

• the rock volume at reservoir condition (BHC) 

• the porosity (<t>) ^ the fluid volume at BHC 

• the initial water saturation (Swi) and net to gross 
(Hu/Ht) ^ the HC volume at BHC 

STOMP = HC volume @ BHC / initial process Bo 




- Process test 


v 1 — o — 

S T 1 P 1 


T = 15°C n 

. ) 


T p P = lbar „ 


O 

L 

T, P- 


Surface 

(separators conditions) 


0 


T 2 p 2 g 

\ 0 


I = Id L 

■ 

P= 1 bar 

1 g 


I — -LD L 

1 


0 — 1‘ 

T, P, P = 1 bar 


o 

T, P 


I = ID L 

1 


T 2 p 2 P = 1 bar 


Reservoir temperature 

1 


Bottom 


0 




r 


L_._, 

P 2 <P 1 


Differential vaporization 


r 


1 bar 


Surface 
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Liquid — Vapor equilibrium 



► Thermodynamics constraints 

• As per thermodynamics, at Liquid-Vapor equilibrium, some of the state functions display 
a minimum, 

• In particular for G = H - T.S (T, P system) ^ dG = - S.dT + V.dP = 0 

• So that, for each component, Fugacity in liquid phase = Fugacity in vapor phase 

► Equation of state (EOS) 

• Since Soave's improvement (1972) the fugacities can be accurately calculated by means 
of an EOS. 

• The most widely used EOS is that of Peng & Robinson (1976): 


RT a(T) 

V-b V(V + b) + b(V-b) 


Fluid studies PVT 
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Peng - Robinson's Equation 



Pure component 



Mixtures 


RT a(T) 

V-b v(v + b)+b(v-b) 


RT an{T) 

V-bm V(V + bm)+bm(V — bm) 


a(T) = a(T c )f(T) 


b = b(T c ) 

a(T c ) = a^ 

' c 


b(T.)-a^ 

r c 

fl L, = 0.4572 


^ = 0.0778 


f(T) 


1+ m(l- Vr”)] 2 


co < 0.49 (nC 10 ) 

m = 0.37464 + 1.54226 co - 0.26992 co 2 


co > 0.49 

m = 0.379642 + 1.48503 co - 0.164423 co 2 + 0.016666 co 3 


b m = Z X I ' b i 


a m 

j=l i=l 

k. = Binary interaction coefficients 
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Fluid model 



► The fluid model is the whole characterization of all the components (Mw, Tc, 
Pc, w, volume shift, Vc, Parachor ...) and the binary interaction coefficient 
matrix (Kij) obtained to perfectly fit the experimental data (PVT report). 

► The fit is obtained by tuning the physical properties of the "heavy" cut(s), the 
Kij (between the lightest and the heaviest components) and the volume shift 
of the lightest and the heaviest component. 

• Tc, Pc, a), Kij to match the equilibrium (saturation pressure, composition), 

• Volume shift to match the volumes, 

• Vc to match the viscosity, 

• Parachor to match the interfacial tension. 


BINARY INTERACTION COEFFICIENTS 

Kij 

GN1 

GN2 

GN3 

GN4 

GN5 

GN6 

GN7 

GN1 


-0,00155 

0,00032 

0,00890 

0,01851 

0,05840 

0,06738 

GN2 

-0,00155 


0,00840 

0,00137 

0,00278 

0,04802 

0,04800 

GN3 

0,00032 

0,00840 


0,00146 

0,00276 

0,02998 

0,02999 

GN4 

0,00890 

0,00137 

0,00146 


0 

0 

0 

GN5 

0,01851 

0,00278 

0,00276 

0 


0 

0 

GN6 

0,05840 

0,04802 

0,02998 

0 

0 


0 

GN7 

0,06738 

0,04800 

0,02999 

0 

0 

0 
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COMPOSITION AND PHYSICAL PROPERTIES 

Name 

ZF % 

Mw 

Tc(K) 

Pc (bar) 

Acent fact 

Vol shift 

Vc (cm 3 ) 

Parachor 

GN1 

71,61 

16,12 

190,31 

45,88 

0,0116 

-0,10 

98,9 

76,6995 

GN2 

11,19 

31,28 

301,25 

49,65 

0,0957 

0 

143,0 

105,4330 

GN3 

7,11 

48,95 

387,44 

40,49 

0,1600 

0 

221,5 

162,1903 

GN4 

3,87 

84,05 

506,62 

32,49 

0,2621 

0 

351,9 

268,0716 

GN5 

2,27 

122,66 

618,15 

26,34 

0,3000 

0 

511,8 

422,4642 

GN6 

3,11 

204,81 

702,45 

18,50 

0,6000 

0 

975,0 

481,6625 

GN7 

0,84 

340,00 

787,25 

17,50 

0,9200 

0,04 

1320,0 

799,5959 
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